Dimensional crossover of the exchange-correlation energy at the semilocal level 
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Commonly used semilocal density functional approximations for the exchange-correlation en- 
ergy fail badly when the true two dimensional limit is approached. We show, using a quasi-two- 
dimensional uniform electron gas in the infinite barrier model, that the semilocal level can correctly 
recover the exchange-correlation energy of the two-dimensional uniform electron gas. We derive new 
exact constraints at the semilocal level for the dimensional crossover of the exchange-correlation en- 
ergy and we propose a method to incorporate them in any exchange-correlation density functional 
approximation. 
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I. INTRODUCTION 

In the Kohn-Sham time-independent density func- 
tional theory (DFT) the noninteracting kinetic energy 
is treated as an exact functional of the occupied or- 
bitalsi and only the exchange-correlation (xc) energy 
Exc has to be approximated. The "Jacob's ladder" 
classification^ of the ground-state density-functional ap- 
proximations for E^c has three complete non-empirical 
rungs: the local-spin-density approximation (LSDA)i, 
the generalized gradient approximation (GGA}^ii, and 
the meta-GGA^. Higher rungs of the ladder require new 
ingredients in order to satisfy more exact constraints^. 
Thus, the meta-GGA has as ingredients the spin densi- 
ties n| and n^, their gradients Vn| and Vn^, and the 
Kohn-Sham (KS) noninteracting kinetic energy densities 
T| and Tj. The local and semilocal density functionals 
(LSDA, GGA, and meta-GGA) give accurate predictions 
of ground-state for atoms, molecules, solids and sur- 
iacei^. They also work for atomic monolayers2ii£ and 
other quasi-two-dimensional (quasi-2D) systems^, but 
they fail badly as the true 2D limit is approachedi^ii^. 
The failure of the semilocal density functionals to de- 
scribe the dimensional crossover of the exact xc func- 
tional can be avoided by using nonlocal modelsi^ii^ such 
as the weighted density approximation^^, or higher rungs 
of the "Jacob's ladder"^. Thus, the fourth-rung hyper- 
GGA, a nonlocal correlation functional compatible with 
exact exchangei^, improves considerably the behaviour 
of semilocal functionals over the whole thickness range 
of the quasi-2D electron gasii. The numerically ex- 
pensive fifth-rung approximations such as the inhomo- 
geneous Singwi-Tosi-Land-Sjolander method (ISTLS)^^ 
and the GW approximation^- are remarkably accurate 
for the description of quasi-2D systemsi^iii. 

The quasi-2D electron gas is experimentally realizable 
in silicon metal-oxide-semiconductor field-effect transis- 
tor (MOSFET) and in the widely used semiconduc- 
tor heterojunctionsi^iSi. Other physical systems with 
strong 2D character are the copper-oxide planes of high- 
temperature superconductors and the electrons bound to 
the surface of liquid heliurr^. 

The paper is organized as follows. In section [TTl we 



present the exact constraints at the semilocal level for the 
dimensional crossover of the exchange-correlation energy. 
In section Hm we construct a simple semilocal functional 
that incorporates these exact conditions, and we test it 
for the quasi-2D uniform electron gas, jellium slabs and 
non-uniformly-scaled hydrogen atom. In section llVi we 
summarize our conclusions. 



II. EXACT CONDITIONS FOR SEMILOCAL 
DENSITY FUNCTIONALS 



A 2D uniform electron gas is described by the 2D 
electron-density parameter r^^ = l/Virn^ 



(Unless otherwise stated, atomic units are used through- 
out, i.e., e'^ = Ti = nie — 1.) Here n^^ is the density 
of electrons per unit area, and kjP represents the mag- 
nitude of the corresponding 2D Fermi wavevector. The 
exchange energy per particle of the 2D uniform electron 
gas i; 
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The correlation energy per particle of a 2D uniform elec- 
tron gas in the high-density limit (r^-° — > 0) is^ 
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A realistic interpolation (which uses Quantum Monte 
Carlo data) between the high- and low-density limits of 
the 2D correlation energy per electron has the following 
formi^i^ 



el^ = 0.5058 



1.3311 
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Similarly to Ref.— , let us consider a quantum well 
of thickness L in the z-direction. In the infinite-barrier 
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model (IBM)^'^ for a quantum well, the KS effective one- 
electron potential is zero inside the well and infinity out- 
side it, such that the KS orbital is 
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for < z < L, 1 > 1, 
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where A is the area of the a;?/-plane, / is the subband 
index, and r| | and k| | are the position and the wavevector 
parallel to the surface. In this model the electrons can 
not leak out of the well, so the true 2D electron-gas limit 
is recovered by shrinking the well. The energy levels of 
this model are^: 



(6) 



When only the lowest level is occupied [Ei f.2D < i?2,o 

which implies L < ■\/3/27rr^-° = £max^^), the density of 
states of this system begins to resemble the density of 
states of a 2D electron gas, the motion in the z-direction 
is frozen out, and the system can be considered quasi- 
two-dimensional. 

By shrinking the z-coordinate without changing the to- 
tal number of electrons per unit area, the system reaches 
the 2D electron gas limit. This process is equivalent with 
a non-uniform scaling in one dimensional, and the 3D 
scaled density isi^ 



{L/X)n{r 



2D\2 



sin^( 



(L/A) 



); 0<z<L/A, 



(7) 

where n\{x,y,z) — \n{x,y,\z) and A is the scaling pa- 
rameter. When A^oo,L/A<< L^a^^, and the 2D limit 
is achieved. The corresponding exchange and correlation 
energies per particle should satisfy the following scaling 
relationsi^ 

lim ^S,K] > -oo ; lim ^E,[nl] > -oo, (8) 

A— ►oo iv A^oo iV 

where iV = n{z)dz. These equations, which start from 
those of Refi^, are not satisfied by LSDA, GGA or meta- 
GGA. 

The GGA exchange-correlation energy per particle of 
our quantum well of thickness L is 

= ( n{z)e^^^{n{z), \/n{z))dz)/N, (9) 

and the meta-GGA exchange-correlation energy per par- 
ticle is 

pMGGA j-L 

i niz)e^i^^^{n{z),Vn{z),T{z))dz) 
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where e^^"^ and e^l.'^'^^ are the GGA and meta-GGA xc 
energies per particle of the 3D system. 



The positive kinetic energy density of the IBM quasi- 
2D electron gas is 
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where = TT{kjPy /{2{L/ \)^) cos^{j^) is the von 

Weizsacker kinetic energy density^. When A ^ oo, t ^ 
t'^ ~ A^ and r-r^ A. Eq. ([TT|) can be well described 
by the Laplacian-level meta-GGA kinetic energy density 
of Ref.26. 

The reduced gradients for exchange { p = 
|Vnp/[4(37r2)2/^n^/^] ~ A^/^) measures the variation of 
the density over a Fermi wavelength and that for corre- 
lation ( t = |Vn|/[4(3/7r)i/Sn^/6] - A^/^) measures the 
variation of the density over the screening length. Both 
tend to infinity when the 2D limit is reached (A oo), 
such that in the quasi-2D electron gas regime, the density 
is rapidly varying almost everywhere. Thus this system 
is not only a challenge for a semilocal density functional, 
but it can give exact constraints (at the semilocal level) 
in the regime where the 3D density and its gradient di- 
verge. 

The exact exchange energy per particle of the 2D uni- 
form electron gas (see Eq. ([1])) can be correctly recov- 
ered by any 3D semilocal density functional approxima- 
tion for exchange energy if in the large gradient limit 
(p oo) the 3D exchange energy per particle behaves as 
a^p-^/^e^'^"^ (see Ref^), where 

= (0.6002V27r)/(3^/2 / dysm''^^(Try)cos-^^^(Try)), 
Jo 

(12) 

was derived in the IBM model using Eqs. ^ and 
Eq. gives ax = 0.5217. The parameter (as well 
as the other results of this section), even if calculated 
using the IBM quasi-2D electron gas, is independent on 
the quantum well potential model along the confinement 
z-direction because the 2D limit (A oo) is not relaying 
on the effective potential model. 

At a meta-GGA level, we can also use the dimension- 
less inhomogeneity parameter— 
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,-unif 



A-2/3, 
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where r"™-'^ = (3/10)(37r2)2/3„5/3 jg ^j^g Thomas-Fermi 
kinetic energy density of the 3D uniform electron gas^-. 
Thus Eq. Il]) can also be exactly satisfied by any meta- 
GGA if in the large gradient limit (p — > oo) the 3D 
exchange energy per particle behaves as bxa^^^e^^^"^ , 
where = 1.947 was found similarly from Eqs. 

dH) and (nni). 

Because r^^ enters in a nonlinear manner in Eq. ([U, 
the GGA level can not describe the correlation energy per 
particle of a 2D uniform electron gas, but it can explain 
the 2D high- and low-density limits. The high-density 
limit of Eq. ^ can be exactly recovered by any 3D 
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GGA that behaves in the large gradient hmit (t oo) 



as 



^GGA 



-0.19 - 0.0497n"^/i2\/t ln(n"^/i2^), (14) 



where t ~ X^^^ is the reduced gradient for correlation. 
The low-density limit of Eq. ([3]) can also be exactly re- 
covered by any 3D GGA that behaves in the large gradi- 
ent limit {t ^ oo) as 

^GGA ^ ^0.40345^5/12^-1/2 + 0.459n5/8f-3/4. (15) 

At the meta-GGA level, Eq.fl]) can be correctly sat- 
isfied by any 3D meta-GGA that behaves in the large 
gradient limit (i — > 00) as 



MGGA 2D 
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where is given by Eq.® and 

r^^ = 0.4173 n-i/3 a'^/^ = 0.6727 r^^ a^^/^. (17) 

Eq. P?)) connects r"^^ with r^^ , showing the impor- 
tance of the a ingredient to the dimensional crossover of 
the exchange-correlation energy, and to the non-uniform 
scaling in one dimension. 



III. CONSTRUCTION AND TESTS OF A 
SIMPLE SEMILOCAL FUNCTIONAL 

The results of Section [IT] can be included in any density 
functional approximation for the exchange-correlation 
energy. For simplicity, we incorporate them in the LSDA. 
Let consider first the exchange part, and define GGA+2D 
and meta-GGA-|-2D as 



^GGA+2D 



LSDA 



{l + /(p)[-l + 0.5217p-i/4]}, (18) 



and 



^MGGA+2D ^ ,LSDA^^ ^ ^ lM7a'^']}, (19) 

where f{p) = 1 for p ^ oo. The simplest approximation 
for f{p) is a step function 



f{p) = lim e{p- x), 

x — >oo 



(20) 



where 9{x) is for x < and 1 for x > 0. This model pre- 
serves all the exact constraints that the local or semilocal 
functional satisfies, and recovers the exchange energy of 
the 2D uniform electron gas in the limit L — (when 
p = oo). However, this approximation does not improve 
the behavior of the semilocal functional in the quasi- 
2D region, and moreover, it gives a discontinuity when 
p — > oo. 

We propose the following simple analytic model for the 
function f{p)^: 



f(p) 



10= + p6 
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where c > is an empirical parameter. Eq. ((2T|) recovers 
the right limit when p oo {f{p —f oo) — 1), and for a 
slowly- varying density, when p is small (p < 1), it behaves 
as f{p) ^ [p^'lO"^ + higher order terms]. This is a good 
feature because it can accurately preserve the behavior 
of the semilocal functional in the slowly-varying limit^. 
When c is large, Eq. dH]) starts to model Eq. pp)) . 

In Fig. [U we show the exchange energy per parti- 
cle of the quasi-2D electron gas with 2D bulk parame- 
ter r^^ = 4 for several density functionals: exact ex- 
change, LSDA, PBE GGA3, GGA-I-2D of Eq. ^ and 
MGGA-h2D of Eq. (HH), using in Eq. ^ three values 
for the parameter c (c =2, 8, and 16). LSDA and PBE di- 
verge when L ^ (A oo). The meta-GGA TPSS^ not 
plotted in Fig. [U has the same behavior as PBE. Both 
GGA-H2D and MGGA+2D perform better than LSDA 
at small thicknesses of the quantum well, and both of 
them recover the exact exchange energy of the 2D uni- 
form electron gas when L = 0. However, we observe that 
GGA+2D with c = 2 (the curve denoted (GGA-f 2D)(a)) 
is the most accurate, and remarcably describes the quasi- 
2D region. When the value of parameter c increases, the 
GGA+2D and meta-GGA-|-2D have the LSDA behav- 
ior over a larger region, and consequently, they are not 
accurate in the quasi-2D region. We also remark that 
GGA-f 2D and meta-GGA-|-2D calculated with c = 8 
(the curves (GGA-f 2D) (b) and (MGGA-|-2D)(b)) give 
a significant improvement over the LSDA in the whole 
quasi-2D region. 

In Fig. [21 we show the exchange energy per par- 
ticle of a thick jellium slab of bulk parameter rf^ — 
2.07. (The bulk parameter defined by the equation n ~ 
3/4:TT{r^^)^ , represents the radius of a sphere that en- 
closes on average one electron). The local and semilocal 
density approximations (LSDA, TPSS, GGA-I-2D, and 
MGGA+2D) show an exponential decay of the exchange 
energy per particle whereas the exact exchange behaves 
as « -l/(4z)32. All the MGGA-I-2D curves have a bump 
in the region where f{p) switches from to 1, and after 
that they are close to the TPSS meta-GGA exchange en- 
ergy per particle. (We recall that TPSS meta-GGA has 
the same large-p behavior as the PBE-GGA) . We observe 
that for jellium slabs (as well as for many 3D systems) 
^GGA+2D > £;LSDA^ whereas i;^^GGA+2i? < £;iSOA, 

thus one can try also a convex combination between 
Eg. p^ and Eg. p^ . However, the construction of an 
accurate 3D and quasi-2D semilocal functional is a dif- 
ficult task^^, and is beyond the purpose of this paper. 

Since p values bigger than 9 are found in the tail of an 
atom or molecule, where the electron density is negligi- 
ble, we can choose c = 8 (such that at p = 9, / = 0.0053 
and for p > 50, / ^ 1). This choice ensures that the 
GGA-f 2D of Eq. ^ and meta-GGA-f 2D of Eq. HH) 
perform similarly with LSDA for 3D systems and make 
a considerable improvement in quasi-2D region, and re- 
cover the exchange energy of the 2D uniform electron gas. 
Thus, from now, all the presented calculations use c = 8 
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FIG. 1: Exchange energy per particle of an IBM quasi- 
2D electron gas of fixed 2D electron density (r^^ — 4), 
as a function of the inverse quantum well thickness 1/L 
(L < Lmax = 15.39). The curves denoted by (GGA-|-2D)(a), 
(GGA-f2D)(b), and (GGA-H2D)(c) are those given by Eq. 
(|18[l. using for /(p) the analytic model of Eq. (|2ip . with 
c = 2, 8 and 16 respectively. The curves denoted by 
(MGGA-h2D)(a), (MGGA-h2D)(b), and (MGGA-H2D)(c) are 
those given by Eq. (I19|l . using for /(p) the analytic model of 
Eq. dHJ, with c = 2, 8 and 16 respectively. While the LSDA 
and PBE diverge, all GGA-I-2D and MGGA-I-2D curves re- 
cover in the limit A ^ cx) the corresponding exchange energy 
(-0.15005) of a 2D electron gas. 



in Eq. dHI). 

Similar with the exchange case, the correlation results 
of SectionlUcan be included in any density functional ap- 
proximation, however, for simplicity we again show them 
for the LSDA case. Thus, we define 



LSDA 



LSDA 



(22) 



where t is the reduced gradient for correlation, f{t) has 
the same form as Eq. (|2ip (with c = 8), and is given 
by Eqs. (g]) and p7)l . 

Figures [3] and |4] show several approximations of the xc 
energy per particle versus the quantum- well thickness L, 
for quasi-2D electron gases of fixed 2D electron-density 



FIG. 2: Exchange energy per particle at position z ver- 
sus z in atomic units, at a jellium slab surface. The 
bulk parameter is rf^ = 2.07, the slab width is d = 
3.2Af and the jellium surface is at 2 = 35.24a.ii. 
(Af ~ 2Tv/k^ is the Fermi wavelength). The calcu- 
lations of the exact exchange, LSDA, TPSS, GGA-f2D 
and MGGA-I-2D use the LSDA Kohn-Sham orbitals. The 
curves (GGA-h2D)(a), (GGA-H2D)(b), (GGA-H2D)(c), and 
(MGGA+2D)(a), (MGGA-f 2D)(b), (MGGA-h2D)(c) has the 
same signification as in Fig. [T] 



parameters: = 4 and — 2/\/3 (as Figs. 1 and 2 of 
Refjii). The ISTLS method^^, a self-consistent approach 
that depends on all occupied and unoccupied KS orbitals, 
is remarkably accurate for any thickness L = L^i^-^/X 
(see Refjii). LSDA and PBE are accurate in the limit 
L — > Lmaxj but they fail badly in the zero thickness limit. 
MGGA+2D and the xc energy per particle of Eqs. (flS]) 
and ^ (e^<^'^+2^ + gMGGA+2D) g^j.^ accurate in the 
limit L Lmax and improve considerably the behavior 
of LSDA when L <^ O.SLmax, approaching the exact 2D 
limit when L ^ 0. 

Let us present our results for the non-uniformly- 
scaled hydrogen atom.^^'^''^, whose density is n\{r) = 
(A/tt) exp(— 2-\/a;2 + + i^z)^)- The exact xc energy 
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FIG. 3: Exchange-correlation energy per particle of an IBM 
quasi-2D electron gas of fixed 2D electron density (r^^ — 4), 
as a function of the quantum well thickness L {L < imax ~ 
15.39). While LSDA and PBE diverge in the 2D limit, 
MGGA+2D and €^°'^+'^° + e'^GA+2D f^^^^ ggj 
(|22|) ') approach the exact 2D limit. 
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FIG. 4: As in Fig. El but now for = 2/^/3 (Lmax = 4.44). 
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(23) 



When A — > oo, this system can model an electron firmly 
bound to a surface. Fig. [5] shows that LSDA and PBE 
fail badly in the extreme oblate case ( A >> 1). The 
meta-GGA TPSS-, not plotted in Fig. [5l has the same 
behavior as PBE. Because in any one-electron system r = 

and a = 0, the MGGA+2D xc energy will approach over, Eq. ^} shows a close connection between r 



FIG. 5: Exchange-correlation energy versus scaling param- 
eter A for the non-uniformly-scaled hydrogen atom. Also 
shown is the exact xc energy (see Eq. (|23p ). 



tionals and keeps as much as possible the 3D accuracy 
of the semilocal functionals. However, in order to obtain 
a good description of the quasi-2D region, the empirical 
parameter c has to be smaller (c w 2) than the optimized 
value (c = 8), but such a choice will significantly mod- 
ify the 3D accuracy of the semilocal functional. In the 
investigation of physical systems with strong 2D charac- 
ter, Eq. pT]) can be seen as an indicator of the quasi-2D 
electron gas regime. Various indicators of the electron 
localization have been constructed for 3D systems (see 
for example Ref i^) , but Eq. P?)) is a better and natural 
choice in the case of quasi-2D uniform gas. Thus, in the 
quasi-2D regions, where r^^(z)a~^/^(z) is constant, we 
can choose c — 2 and in the other regions we can choose 
c = 8 (or even c — oo). Such a parametrization of c 
can be more useful in applications than the use of our 
optimized value for c (c = 8). 



IV. CONCLUSIONS 

In summary, we have shown that the dimensional 
crossover (from 3D to 2D) of the exact xc energy can be 
significantly improved at a meta-GGA level, and we de- 
rive new exact constraints (see Section |lT]) using an IBM 
quasi-2D electron gas. Same results can be obtained us- 
ing the parabolic quantum well of Ref^i^ because the 2D 
limit is independent on the quasi-2D electron gas model. 
Thus, a 3D meta-GGA that requires input from the 3D 
uniform electron gas-''', can describe a 2D system using 
only the highly nonlocal region where |Vn| oo. More- 

2D 



slowly zero in the limit A — > oo. The xc energy of Eqs. 
(Ull) and ^ (^E^GA+2D ^^MGGA+2D-^ ^^^^ accurate 

at large values of A. 

In this section we have proposed a method that incor- 
porates the 2D limit of the electron gas in semilocal func- 



and the positive noninteracting kinetic energy density r 
in the case of the quasi-2D electron gas. 

We propose a simple approach to incorporate the di- 
mensional crossover constraints in any local or semilocal 
density functional approximation for the xc energy, and 
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we present it in the case of LSDA. However, future work 
is needed to construct an accurate meta-GGA that sat- 
isfies the dimensional crossover constraints. 

The non-uniform scaling in one dimension is closely 
related to the quasi-2D electron gas. The non- uniformly 
scaled hydrogen atom in the oblate case (A > 1), an im- 
portant and hard test for the density functionals as well 
as a model for an electron bound to a surface, can be well 



described by our simple modified LSD As. Thus we hope 
that this work can be useful not only for investigation of 
physical systems with strong 2D character, but also in 
developing more accurate density functionals. 
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sor J.M. Pitarke for many valuable discussions and sug- 
gestions. L.A.C. acknowledges NSF support (Grant No. 
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